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1 Introduction 

The SQUID magnetometer system with ultra-high 
magnetic sensitivity has made it possible to measure 
the very weak magnetic field generated by the 
activity of the brain and the heart. Biomagnetic 
measurement of the magnetic field perpendicular to 
the body surface is widely used. We had developed 
a 39-channel SQUID system consisting of thirteen 
3-D second-order gradiometers (12 measurement 
positions (4x3) and a noise measurement) for vector 
measurement of the magnetoencephalogram (MEG) 
and the magnetocardiogram (MCG). This SQUID 
magnetometer can detect magnetic field components 
perpendicular to the body and tangential to the body 
simultaneously [1], We carried out 3-D vector 
measurements of MEG with a mixed auditory 
evoked field (AEF) and somatosensory evoked field 
(SEF) overlapping in time, and MCG with normal 
subjects and an IRBBB (incomplete right bundle 
branch block) patient were also measured [2], In 
order to discuss the frequency characteristics with 
time course, a short-term spectrum of biomagnetic 
signals were calculated by an autoregressive (AR) 
model, and the power spectrum was analyzed (time- 
frequency analysis). 

2 Methods 

2.1 Biomagnetic measurement 

We carried out 3-D measurement of MEG and MCG. 
Fig. 1 a) showed a 3-D second order gradiometer, 
which was orthogonal wound with NbTi-Cu wire on 
a rectangular solid 3x3x6 cm. 

1. MEG measurement was the mixed AEF and 
SEF overlapping in time. AEF was elicited by 1 
kHz tone bursts of 20 ms duration to the right 
ear. SEF was elicited by electric pulses of 0.2 
ms duration with 5 mA to the median nerve of 
the right wrist. In the mixed AEF and SEF, 
stimulation to the median nerve was delivered 
80 ms later than to the right ear. Fig. 1 b) 
showed coordinate system and measurement 
positions for MEG. There were forty-eight 


measurement positions (4 scanning, fixed circle 
in Fig. 1 b)) on the scalp and a band-pass filter 
was used in the range of 0.1 to 40 Hz. All 
magnetic data were averaged for 300 

measurements at each position. 

2. MCG measurements were carried out to normal 
subjects and an IRBBB patient, and ECG (lead 
II) was measured simultaneously. Fig. 1 c) 
showed coordinate system and measurement 
positions for MCG. There were forty-eight 
measurement positions (4 scanning, fixed circle 
in Fig. 1 c)) on the chest wall and a band-pass 
filter was used in the range of 0.1 to 80 Hz. All 
magnetic data were averaged for 50 

measurements at each position. The reference 
signal for averaging was R-wave of ECG. 

2.2 Time-Frequency analysis 

Time-frequency analysis of the biomagnetic data 
was calculated by the AR model [3],[4], This AR 
model is effective in respect of frequency resolution 
in short-time spectrum. The analyzing time was 
50ms and/or 100ms. The order of AR model was 
decided by final prediction error criterion (FPE) 
between 15 and 20. The short-time spectrum is 
calculated by equation (1). 


1 + Y_i a k exp(-/2;z/M/)| 

Where B(f)\ power spectrum of magnetic field, <7;: 
variance of white noise, af : coefficient of the AR 
model, p: order of AR model, At : sampling 
interval. 

3 Results 

3.1 MEG measurement 

Fig. 2 shows the waveform (Fig. 2 a)) of the mixed 
AEF and SEF overlapping in time measured at 
position F2 (Br component), the time-frequency 
analysis (Fig. 2 b)), and the intensity distributions of 





Figure 1: a) A 3-D second order gradiometer. b) 
Coordinate system and measurement positions for 
MEG. c) Coordinate system and measurement 
positions for MCG. 

power spectrum. Fig. 2 c)-e) and f)-h) show the 
intensity distributions of 15Hz power spectrum at 
latencies of 38 ms (analysis time: 38-88 ms) and 150 
ms (analysis time: 150-200 ms) respectively. These 
figures show the intensity distribution of power 
spectrum obtained from each measurement position. 
In Fig. 2 c)-e), it can be seen the response at G4 by 


auditory stimulation before median nerve is 
stimulated. B0 component of Fig. 2 d) shows the 
peak at G4. The peak of tangential components (B0, 
B(])) corresponds to the source location [5], Namely, 
the source can be estimated at G4. The position G4 
corresponds to the location of auditory area at the 
left hemisphere [6], In Fig. 2 f)-h), it can be seen the 
response at D3 by median nerve stimulation. The 
position D3 corresponds to the location of 
somatosensory area at the left hemisphere [7]. 

Fig. 3 also shows the waveform of the mixed AEF 
and SEF overlapping in time measured at position 
F2 (Br component), the time-frequency analysis, and 
the intensity distributions of distributions of 35 Hz 
and 15 Hz power spectrum at latency of 50 ms 
(analysis time: 50-100 ms) respectively. B0 
component of Fig. 3 d) shows the peak at D3. 
Therefore, the source of 35 Hz component can be 
estimated at D3. In Fig. 3 c)-e), it can be seen the 
response by median nerve stimulation, by referring 
35 Hz component shown in Fig. 2 c)-e). B0 
component of Fig. 3 g) and B(|) component of Fig. 3 
h) show the peak of 15 Hz component at G4 
receptivity. In Fig. 3 f)-h), these can be seen the 
response by auditory stimulation, by referring 15 Hz 
component shown in Fig. 2 c)-e). 

By considering time-frequency analysis, the 
separation of the multiple sources is possible. These 
results show that time-frequency analysis is a useful 
method for estimating multiple sources. 



Figure 2: a) The waveform of the mixed AEF and SEF overlapping in time measured at position F2 (Br 
component), b) The time-frequency analysis. c)-d) The intensity distributions of power spectrum of 15 Hz 
at latency of 38 ms ( c):Br, d):B9, e):B(j>). f)-h) The intensity distributions of power spectrum of 15 Hz at 
latency of 150 ms (j):Br, g):B0, h):B(j)). 

















Figure 3: a) The waveform of the mixed AEF and SEF overlapping in time measured at position F2 (Br 
component), b) The time-frequency analysis. c)-d) The intensity distributions of power spectrum of 35 Hz 
at latency of 50 ms ( c):Br, d):BO, e):B<j)). f)-h) The intensity distributions of power spectrum of 15 Hz at 
latency of 50 ms (f):Br, g):BO, h):B<5). 


3.2 MCG measurements 

Fig. 4 shows the waveform (Fig. 4 a), c), e)) of the 
normal subject and time-frequency analysis (Fig.4 b), 
d), f)). In Fig. 4 b) Bz component, it is shown that 
the component of about 14 Hz dominates for 450 ms. 
It is also shown that the component of about 20 Hz 
dominates for 460 to 490 ms and about 15 Hz 
dominates for 510 ms. Fig. 4 d) Bx component 
shows the same pattern as Bz component. In Fig. 4 
f) By component, it is shown that the component of 
about 17 Hz dominates for 450 ms and about 24 Hz 
dominates for 460 to 480 ms. These time courses 
show the ventricular depolarization (the QRS 
complex). 

Fig. 5 shows the waveform (Fig. 5 a), c), e)) of the 
IRBBB patient and time-frequency analysis (Fig.5 
b), d), f)). In Fig. 5 b) Bz component, it is shown 
that the component of about 27 Hz dominates for 
470 ms, and about 15 Hz dominates for 490 to 510 
ms, and it differs from normal subjects. Especially, 
the component of about 10 Hz for 550 ms (S-wave) 
appears only in the IRBBB patient. Fig. 5 d) Bx 
component and f) By component show the same 
pattern as Bz component.. 

Time-frequency analysis can show the change of 
frequency component on the MCG waveform with 
the time. The frequency of R wave of the normal 
subjects and the IRBBB patient were about 20 Hz 
and 27 Hz respectively. And each frequency was 
different. In the IRBBB patient (Fig. 5), it is shown 


that the component of about 10 Hz dominates for 
550 ms. 

4 Conclusion 

In the MEG, the mixed AEF and SEF overlapping in 
time, spatial distributions of power spectrum of 
15Hz and 35Hz at 50-100m s clearly showed that 
two sources existed, by referring tangential 
component. The latencies corresponded to two 
sources were N100 of AEF and N20 of SEF. The 
tangential component of 35Hz indicated an extreme 
at the position corresponded to somatosensory area. 
The tangential components of 15Hz indicated an 
extreme at the position corresponded to auditory 
area. Two sources were separated with different 
frequency by using time-frequency analysis. 

In the MCG, time-frequency analysis to an IRBBB 
patient showed the extreme at S-wave with about 10 
Hz. However, this frequency component didn’t 
show in normal subjects. It is found that time- 
frequency analysis is a useful method for estimating 
multiple sources and understanding the fine 
structure of biomagnetic signals. 
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